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Abstract 



In strong perpendicular magnetic fields double-quantum-well systems can 
sometimes occur in unusual broken symmetry states which have interwell 
phase coherence in the absence of interwell hopping. When hopping is present 
in such systems and the magnetic field is tilted away from the normal to the 
quantum well planes, a related soliton-lattice state can occur which has kinks 
in the dependence of the relative phase between electrons in opposite layers on 
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the coordinate perpendicular to the in-plane component of the magnetic field. 
In this article we evaluate the collective modes of this soliton-lattice state in 
the generalized random-phase aproximation. We find that, in addition to the 
Goldstone modes associated with the broken translational symmetry of the 
soliton-lattice state, higher energy collective modes occur which are closely 
related to the Goldstone modes present in the spontaneously phase-coherent 
state. We study the evolution of these collective modes as a function of the 
strength of the in-plane magnetic field and comment on the possibility of 
using the in-plane field to generate a finite wave probe of the spontaneously 
phase-coherent state. 
Pacs: 73.20.Dx,64.70.Rh,73.20.Mf 

I. INTRODUCTION 

Recently, there have been a number of experimental 0] and theoretical |0| investigations 
of the integer and fractional quantum Hall effects in double-quantum- well systems (DQWS). 
The additional degree of freedom due to the presence of the second well makes the phase 
diagram of the two-dimensional electron gas in the DQWS very rich. Technological progress 
has opened the possibility of producing DQWS's of very high mobility where the spacing 
between the wells is small enough (rf ~ 100 A ) to be comparable with the typical elec- 
tron spacing within a layer. In this regime a DQWS may differ qualitatively from a pair 
of isolated single-layer systems. In particular, theoretical work based on several apparently 
different points of view has established 0-§l that for sufficiently small separations, a novel 
broken symmetry can occur in DQWS's in a strong perpendicular magnetic field in which 
phase coherence exists between electrons in different layers even when no tunneling occurs. 
This broken symmetry is favored by interlayer electron-electron interactions because good 
interlayer correlations are guaranteed by the Pauli exclusion principle when the phase differ- 
ence between the two layers is fixed. (If the phase relationship is fixed, two electrons cannot 



share the same planar coordinate even if they are in opposite layers.) The broken symmetry 
occurs at strong perpendicular magnetic fields because Landau level degeneracy allows the 
symmetry to be broken without a kinetic energy cost. The broken symmetry is expected 
to exist in the ground state for layer separation, d, smaller than a critical layer separation, 
dc- Moreover, for d < dc the system is expected |^J^ to have a finite temperature phase 
transition of the Kosterlitz-Thouless type. 

For realistic DQWS's with d < dc there is always a finite amplitude for tunneling be- 
tween the two layers of the DQWS. Tunneling favors a particular phase relationship between 
electrons in opposite layers so that the symmetry discussed above is explicitly broken by the 
tunneling term in the Hamiltonian and the Kosterlitz-Thouless phase transition associated 
with its breaking does not occur. However, recent experiments [0 have demonstrated that a 
second related quantum phase transition occurs in DQWS's when the magnetic field is tilted 
away from the normal to the layers of the DQWS. When the component of the magnetic 
field which is parallel to the electron layers is finite, electrons must acquire an Aharonov- 
Bohm phase when they move around paths which enclose fiux from the parallel field. In one 
convenient gauge, the Aharonov-Bohm phase appears through phase factors in the inter- 
layer tunneling amplitudes which depend on the planar coordinate which is perpendicular 
to the in-plane component of the magnetic field. The tunneling term in the Hamiltonian 
then favors a phase relationship between electrons in opposite layers which follows the phase 
of the tunneling amplitude. When the phase relationship is spatially dependent it becomes 
possible for electrons in opposite layers to share the same planar coordinate so that the 
interaction energy increases. 

According to Yang et al. [HH,^ the phase transition seen by Murphy et al. results 
from this competition between tunneling energy and interaction energy in a parallel field. 
In the simplest view, the state on the small parallel field side of the transition is one where 
the phase relationship between the layers follows the phase of the tunneling amplitude 
while the state on the large parallel field side of the transition is one in which the phase 
is independent of the spatial coordinate and the symmetry-breaking tunneling term has 



no effect on the ground-state energy of the DQWS. A more accurate picture of this state 
and of the phase transition follows from the effective Hamiltonian derived by Yang et al. 
in which the local interaction energy cost of spatial dependence in the phase relationship 
between the layers is approximated by a gradient term in the energy density. Minimizing the 
energy in a parallel field then gives rise to a sine-Gordon equation. The state on the large 
parallel field side of this transition is a soliton-lattice state H,^ which has its translational 
symmetry broken by the introduction of periodic kinks in the phase relationship between 
electrons in the two wells of the DQWS as the parallel magnetic field is increased. The 
effective Hamiltonian derived by Yang et al. may be mapped to thoroughly studied models 



(see for example Refs. |jTO|JTT|) of commensurate-incommensurate phase transitions with the 
parameter which tunes the transition proportional to the parallel component of the magnetic 
field (see below). Motivated by this mapping, we refer to the state on the weak parallel field 
side of the transition as the commensurate state and to the state in the absence of tunneling 
as the incommensurate state. 

In this paper we use a generalized random-phase approximation (GRPA) to study the 
response functions and collective modes of the sohton-lattice state. As discussed elsewhere ^ 
and in further details below, the GRPA is expected to be qualitatively correct for DQWS's 
at total Landau level filling factor u = 1. It is at this filling factor that the physical 
effects associated with spontaneous phase coherence are most easily observed. The use 
of a microscopic theory frees us from the gradient approximation and, more importantly, 
allows us to study higher energy excitations which may not be accurately rendered by the 
effective Hamiltonian. In Section II, we detail the microscopic Hamiltonian on which we 
base our calculations. We appeal to the strong perpendicular magnetic field in limiting the 
Hilbert space to the lowest Landau level in each layer and assume that the electrons are 
completely spin polarized. GRPA response functions are readily calculated from the Hartree- 
Fock approximation ground state of the DQWS. In Section HI, we outline the Hartree- 
Fock formalism which allows us to describe both the simpler commensurate state and the 
more intricate soliton-lattice state. Numerical Hartree-Fock approximation results for the 



ground state are compared with expectations based on the gradient approximation effective 
Hamiltonian and with the theory of commensurate-incommensurate phase transitions. At 
strong parallel magnetic fields corrections are quantitatively important. We discuss results 
for the period of the soliton-lattice, the shape of individual solitons, and the ground-state 
energy. We show explicitly that the soliton-lattice state always has lower energy than the 
incommensurate phase, although the energy difference is extremely small at very strong 
parallel magnetic fields. We also calculate the critical strength of parallel magnetic field 
at which the phase transition between commensurate and soliton-lattice states occurs. In 
Section IV we describe the GRPA calculations. The response functions of the soliton-lattice 
state are complicated by its broken translational symmetry but we are able to derive formal 
results which are convenient for numerical calculations by taking advantage of simplifications 
which are due to the truncation of the Hilbert space to the lowest Landau level. The 
formalism is a generalization to double-layers of one developed [Q to describe the Wigner 



crystal state of single layer systems in strong perpendicular magnetic fields and has been 
applied previously |TB| to double-layer systems in the case where there is no parallel field. We 
find that the lowest energy collective modes are Goldstone modes associated with the broken 
translational symmetry of the soliton-lattice state. Collective modes occur at higher energies 
which are closely related to the Goldstone modes of the phase-coherent broken-symmetry 
state in the absence of tunneling between the layers. We comment on the possibility of 
using a parallel magnetic field to generate a finite wave vector probe of the coherent broken 
symmetry state. Section V contains a brief summary of our results. 

II. MODEL HAMILTONIAN FOR DQWS'S IN A TILTED MAGNETIC FIELD 

We consider a symmetric DQWS subjected to a tilted magnetic field. We are interested 
in studying the ground-state properties of this system as a function of both the parallel 
component of the magnetic field and the separation d (from center to center) between the 
wells at filling factor z/ = 1. We assume throughout this paper that the perpendicular 



component of the magnetic field is kept constant. The total field is taken as B = i?|| y + B±7, 
and, in the Landau gauge, the vector potential A = (0, B^ x, —B\\x). We consider only the 
case of an unbiased DQWS: the subband energies in the two layers are identical and, in the 
ground state, the charge is equally distributed between the two layers. 

Taking the perpendicular component of the magnetic field to be strong, we restrict the 
Hilbert space to the first Landau level and assume the electron gas to be fully spin polarized. 
The many-body Hilbert space then consists of occupation number states based on a set of 
single-particle orbitals labeled by a layer index i and a guiding center index X: 
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where i'^ = hc/eBj_ is the square of the magnetic length for the perpendicular component 
of the magnetic field and Xi{^) with i = R,L is the envelope wave function of the lowest- 
energy state centered on the right or left well. We can define a separate magnetic length 
associated with the parallel component of the magnetic field, i?, = hc/eB\\. In the absence 
of interactions and tunneling between the wells these single particle orbitals are eigenstates 
of the single-particle Hamiltonian. 

For the sake of limiting the number of parameters characterizing our DQWS, we make 
a narrow well approximation i.e. we assume that the width, 6, of the wells is small (6 << 



d) and treat interlayer hopping in a tight-binding |T^ approximation. The single-particle 
problem is then characterized by the interlayer distance d and a hopping integral t: 

Ho = T.^ ixc^,x - t-E {e-''''^'' 4,xCL,x + h.c.) . (2) 

i,X X 

In this approximation, the parallel component of the magnetic field is responsible for a 
reduction [T^ in the magnitude of the hopping integral (t = te~ n = te~'^ ^ Z^) and 



for the introduction of a guiding-center-dependent phase factor, e~* ' ii = e"**^^. e is the 
identical subband energy of the quantum wells (including the lowest Landau-level energy 
huJc/2, where Uc = eB±/mc) and Q is defined by 

Q = ^oci?||. (3) 



The parallel component of the magnetic field enters only in the non-interacting part of the 
Hamiltonian but, as we will see, the competition it produces with the interaction energy 
leads to interesting many-body physics. 

III. HARTREE-FOCK GROUND STATES IN A PARALLEL FIELD 

A. Hamiltonian of the interacting 2D electron gas 

In order to evaluate the GRPA response functions we must first complete self-consistent 
Hartree-Fock calculations for the ground state of the DQWS. The order parameters that 
describe the commensurate and soliton-lattice states are defined below and we compute 



them using the formalism developed in Refs. |13] and |12]. Although this approach is not 



the simplest approach one could use in the present context, it has the advantage of being 
easy to implement numerically. More importantly, the results we derive here will lead to a 
remarkable simplification of the GRPA response function calculations. 
We start be defining the operator 

P^M) = - E e-^"^^-"^'^^^^/^ 4x9,x+,„.^, (4) 

9 X 

where g = S/2'kI'^ is the degeneracy of the Landau levels [S is the area of the system). 

The single-particle Hartree-Fock Hamiltonian can be expressed in terms of this operator. 

Defining Q = Qx we obtain 

Hhf = ^ep(O) - gt{paL{Q,) + Plr{-Q,)) 

(5) 
+ 9 (^) E*,i Eq [HiijM) < Pn(-q) > Pii(q) - Xiijjiq) < Pij{-q) > pji{q)] , 



where 



Hu^^icl) = K(q) = (^) e-^'^y^ 



are respectively the Hartree intrawell and interwell interactions and the Fock intrawell 
and interwell interactions. These functions of wave vector depend only on the interlayer 
separation and are plotted in Ref. [|1^]. Jo is a Bessel function of the first kind and 
p(0) = Prr{0) + Pll{0) is the "density" operator. In deriving these results we have ne- 
glected processes in which electrons are scattered from one well to the other as the result 
of their mutual interactions. This assumption is consistent with the tight-binding approx- 
imation. In realistic systems, small terms of this type will be present but they do not add 
anything new to the physics. 

These equations are very general and allow both positional and coherence broken sym- 
metries. The set of quantities < pij (q) > describes all possible states of the DQWS and we 
use them as order parameters: < PRR{q) > and < Pll{^) > describe density fiuctuations 
in each well while the interwell coherence is refiected in a non-zero value of the interwell 
order parameters < PRL{<i) > and < PL_R(q) > • The states that we consider here are ones 
where the electron density is uniform in each well, but the interwell order parameters are 
allowed to vary in space along the x axis. We do not allow variation in the interwell order 
parameters along the y direction. In our notation, we restrict ourselves to solutions of the 
Hartree-Fock equations where 

< pRRiO) > = < Pll(O) >= i, 

< pRLiqx, 0) > , < pLRiqx, 0) >^ 0, 
and all other order parameters are zero. It seems clear that the lowest-energy solution of 

the Hartree-Fock equations satisfies these restrictions for the situations of physical interest 

to us. The restriction that < Prr{0) >=< Pll{0) >= 1/2 is enforced by the large Hartree 

energy cost of having unequal densities in the two layers. The only spatial dependence 

which we allow is that driven by the hopping term in the Hamiltonian which depends on 

the X coordinate. In the absence of a parallel field, we know [|13| that the minimum-energy 

solution of the Hartree-Fock equations is spatially uniform at small layer separations but 

breaks translational symmetry by forming a Wigner crystal state at sufficiently large layer 

separations. Once the Wigner crystal state is formed, interlayer coherence is quickly lost. 
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For u = 1 the occurrence of a Wigner crystal state is believed to be an artifact of the 
Hartree-Fock approximation. In any event, our present interest is in the regime of small 
layer separations where we do not expect spatial dependence other than that driven by the 
parallel field; here we are confident that the lowest-energy Hartree-Fock state satisfies the 
restrictions of Eq.( |^. With these restrictions the Hartree-Fock Hamiltonian simplifies: 

Hhf = g{e-^, {^) H(0)) p(0) - gtipRLiQ) + Pl/?(-Q)) 

\°) 
- 9 {^) Eg, Vd{qx) [< pRii-qx) > pLR{qx)+ < pLR{-qx) > pRhiqx)] , 



where V{,(0) = wf . To arrive at Eq.(Q), we have added a neutralizing positive background; 
the energy e appearing in Eq.(H) now depends on the location of this background in three 
dimensional space. For constant perpendicular magnetic field, the first term in the right- 
hand-side of Eq.@ is a constant which we absorb into the zero of energy in what follows. 

B. Equation of motion of the single-particle Green's function 

The Hartree-Fock Hamiltonian depends on order parameters which must be determined 
self-consistently. We now derive explicit self-consistent equations by considering the equa- 
tion of motion of the single-particle Green's function defined by G'jj(X, X,r) = — < 
Tcj^x('7")c j^'(O) > . Within the restricted Hartree-Fock approximation discussed above, 
the single-particle Green's function is diagonal in its guiding center indices and we can work 



15 with a one-component Fourier transform defined by 



G^A<i-, r) = -Y. Gr^X. X, r)e-^'^-^. (9) 

9 X 

The order parameters are related to this Green's function by < Pij{qx) >= GjiiqxjT = 0^). 

The equation of motion of the Green's function is derived using h-^(- • ■) = [H — pN, ■ ■ •] 

where p is the chemical potential of the electron gas and we set ^ = 1. We find that 

{iuJn + p)GRR{qr„ Un) + tGlRiqx + Q, tUn) = ^fe - Eg' U*{q'x)GLR{qx + q'cc, Wn) , , 

(10) 
{iUn + p)GLR{qx,(^n)+tGRR{q^-Q,uJn) = - Y.q'^U{q'^)GRR{q^ - q'^, u^), 



where the effective Hartree-Foclc mean field is due to the interlayer exchange energy and is 
given by 

U{qx) = (t^J MQx) < pRL{qx) > . (11) 

Our restricted Hartree-Fock equations do not couple different guiding centers. Since we 
require that the charge density in each layer be uniform and independent of position, it 
follows that for the Hartree Fock eigenfunctions at guiding center X, the magnitudes of the 
coefficients multiplying the basis functions localized in left and right wells must be equal. 
The only freedom of physical relevance is the relative phase, 0{X), multiplying the basis 
functions. In terms of this quantity, 

< PRLiq.) >= -E < P(^) > ^-''^'' = -E 11^''^''^] e-^"^"", (12) 

9 X 9 X ^^ ^ 

where < p{X) >=< c\ix'^l,x > is the interwell order parameter in X space. 
At this point, it is more convenient p|,p to define a new phase 0{X) by 

e{X) = e{X) - QX. (13) 

In the commensurate state, 0{X) = 0, while in the soliton-lattice state it is a periodic 
function of X. It is advantageous to make use of this periodicity and redefine the interwell 
order parameter as 

< pnLiq.) >=< PRLiq. + Q) >= 7^ E e-^'^^^e*^(^), (14) 

^fi' X 

which is the Fourier transform of < p{X) >=< p{X) > e~^'^^ = |e*^*^"^^ . 

In our notation the commensurate-incommensurate phase transition of Ref. P] is de- 
scribed in the following way. The commensurate phase (C phase) at weak parallel magnetic 
fields has 0{X) = QX; the phase of the wavefunction follows the phase of the hopping 
matrix element just as it would for non-interacting electrons. For this state 0{X) = and 
so only < Prl{0) >= | is different from zero. The commensurate phase has an energy per 
particle 
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which increases monotonically with magnetic field. In the opposite hmit of strong parallel 
magnetic fields, the phase does not follow the imposed periodicity and we expect 6{X) = 0. 
Thus, 0{X) = —QX or, equivalently, only < Prl(— Q) >= ^ is non-zero. The energy per 
particle in the incommensurate phase (I phase) is 

and is clearly independent of the parallel magnetic field. There is a critical parallel magnetic 
field B\\^(^c^i), at which the energy of the commensurate state rises above the energy of the 
incommensurate state. As pointed out in Ref. P], however, before this critical magnetic field 
is reached it becomes energetically favorable to introduce kinks (solitons) where 0{X) changes 
by 2tt] in the ground state the solitons are arranged periodically. We refer to this state as the 
soliton-lattice state (S phase). For very strong parallel magnetic fields the solitons are smooth 
and spaced by 27i/Q so that 9{X) ^ —QX and the soliton-lattice state asymptotically 
approaches the incommensurate state. In Ref. |^ this phase transition is described using 
a gradient approximation for the exchange energy which allows the ground-state problem 
to be mapped to well-studied models of commensurate-incommensurate phase transitions 
10| , p!Tl| . In the Appendix, we briefiy review the most relevant results from these models and 



discuss the relationship between the gradient-approximation theory for the ground state 
and the Hartree-Fock approximation. At the critical parallel field for the transition between 
commensurate and soliton-lattice states, B\\^(^c'^s), the system admits a single soliton and 
0{X) is a solution of the sine-Gordon equation. (See the appendix for further details.) As 
the parallel magnetic field increases above this transition, the solitons proliferate and it is 
necessary to take their interactions into account. This is best done numerically p!0| , pJ]| even 
in the gradient approximation. Below we describe our numerical Hartree-Fock calculations 
for the soliton-lattice ground state. 

In the ground state of the soliton-lattice phase, we have the periodicity e*^^"^"'"^^) = ^t^ix) 
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where Ls, the period of the sohton-lattice, is chosen to minimize the energy. It follows that 
if we define a soliton wave vector by 

271 



Qt 



(17) 



the only non-zero order parameters are < PblI'^Qs) >, n = 0, ±1, ±2, . . .. (We can choose 
the origin for the x coordinate so that these order parameters are real.) We can simplify 
the notation of Eqs.(0) by redefining the interwell single-particle Green's function as 



GlRiqx, UJn) = GLR^qx + Q, UJn) 



(18) 



Defining GLR{n,uJn) = GLRinQs,uJn),GRR{n,uJn) = GRR{nQs,uJn), u{n) = Sn,o and the 
real matrix 



F{n, m) = t6n,m + ( ~« 1 ^d{{n - m)Qs + Q) < pRhin - m) >, 
we have, in an obvious matrix notation, 

{iUn + p)GRR{Un) + F'^GlrM = U, 
{iUn + fi)GLR{uJn) + FGRR{uJn) = 0. 

These equations are readily decoupled: 



[iun + At) I - B 



GRR^iOn) 



{iuJn + yU) M, 



{iuJn + p) I - B\ GiRiuJn) = ^Fu, 

where the real and symmetric matrix B is defined by 

B = F^F = FF^. 



(19) 



(20) 



(21) 



(22) 



B can be diagonalized by an orthogonal transformation: B = UQU'^ where Q is the 
diagonal matrix of the eigenvalues of B with elements Qi. Solving for the order parameters, 
we find at finite temperature 



(PRRin)) = 2 E Un.Uj;, [f{E,,.) + f{E,,+)] , 



(23) 
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and 

(m(n)) = \ E,. '"""fe;^''"^ [/(i?.-) - /(i?.+)] , (24) 

where f{x) = ^^^^i^) is the Fermi function and 

E,,± = ±^. (25) 

In order to avoid density fluctuations at T = OK, i.e. in order to have (pii(O)) = 
i, (pM(n 7^ 0)) = (i = -R, L) , we must have fi{T = 0) = so that the Fermi factors 
equal one for the negative energy solutions and zero for the positive energy solutions. Then 
{pii{n)) = ^J2j UnjUjQ = ^6nfl and, in matrix notation, 

(Prl) = ^UQ-^^U^Fu, (26) 

where f2^'^ represents the diagonal matrix of fi/ values. From this last equation, we can 
directly show that the interwell order parameter obeys the sum rule 

En\{pRL{n))\^ = {pRhf {Prl) , 

= \xy^F^B-^Fn, 

_ 1 

4' 

Eqs. (p2[) , (p5[) , (|26|) represent a set of self-consistent, coupled equations that we must solve 
in an iterative way until a stable solution is found for a given value of Qs, Q and t. The 
calculation must then be repeated for different values of Qs until we find the solution which 
minimizes the HF ground-state energy per particle 

E = -2t< puLiO) > - (^] E M^Qs + Q) {pRL{n)f . (28) 

Given the order parameters, we can easily find 6{X) from the equation 

< Prl{X) >= ^e'^(^) = Y: e^"^^"" < PRL{n) > • (29) 

The Hartree-Fock eigenvalues E± {X) (see appendix) are related to the order parameters by 
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E±{X) = ± 



t + E VciinQs + Q)< PRL{n) > e^'^sx 



(30) 



Commensurate and incommensurate states also represent extrema of the Hartree-Fock 
energy functional and therefore are self-consistent solutions of the Hartree-Fock equations. 
For commensurate and incommensurate states respectively we find that 

E^iX) = ± (i+ Iv.iQ)^ , (31) 



and 



Ef{X) = ± 



t + ^Kz(0)e-*«^ 



(32) 



C. Numerical results for the Hartree-Fock Approximation Ground State 

In Fig. 1, we plot the energies of the C, I and S phases for (a) t/ {(? jt) = 0.01, d/i = 1.0 
and (b) t/ (e^/£) = 0.005, d/i = 1.877. The dotted line represents the gradient approxima- 
tion for the commensurate state energy 

Ecigrad) ^ -t - \ i^\ V^ + Tl^ psQ^ (33) 

The stiffness ps is easily obtained from Vd{Q) using Eq.( |A14| ). We find for (a) vrps = 
0.0195 (e^/£) and (b) ixps = 0.00751 [e'^/i). The gradient approximation is very close to the 
HF energy for Eq- Analytic results from the gradient approximation for the the values of Q 
at which the commensurate and incommensurate states are equal in energy {Qc^i) and at 
which the soliton state first becomes stable {Qc^s) are also in good agreement with the the 
HF values. Numerically, we find that for (a) Qc-*s/Qc^i = 0.91, and for (b) Qc^s/Qc^i = 
0.95 while in the gradient approximation, this ratio is Qc-^s/Qc^i = 2v^/7r = 0.90. 

We also see in Fig. 1 that the energy of the soliton phase is always lower than the energy 
of the incommensurate phase so that this phase is never thermodynamically stable in the 
Hartree-Fock approximation. Of course, as Q becomes very large, the order parameters in 
the soliton phase become asymptotically close to those of the / phase and the energies of the 
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two phases becomes nearly identical. One should also notice that the difference in energy 
between the soliton and incommensurate or commensurate phase is of the order of 1% at 
most! 

The behavior of 0{X) over one period of the soliton lattice is shown in Fig. 2. As 
expected, the width of the soliton is very large right at the transition and then decreases 
as the parallel magnetic field increases. In the large Q limit, 0{X) -^ —QX + vr, which is 
the incommensurate solution (the constant vr here comes from the fact that we have made 
an overall phase convention which forces the soliton solution to vary from 27r to over the 
interval {-Ls/2,Ls/2)). 

The dotted line in (a) and (b) represents the sine-Gordon solution 0{X) = 
4tan"^ \e~^/^ 



, where ^ = J27ipsi'^/t (see Eq.( |A19| )). Fig. 3, is a plot of the soliton wave 



vector Qs = 27^/ Ls as a function of Qi. As the parallel magnetic field increases, the soliton 
wave vector asymptotically approaches Q as required to approach the incommensurate state 
solution. 

Fig. 4 shows the shape of the upper band Ej^{X) over one period Ls for two extreme 
values of QC. The dispersion of the soliton band structure seen in (a) weakens as the parallel 
magnetic field is increased (see (b)). In the Hartree-Fock approximation the charge gap at 
z/ = 1, the minimum energy to make a particle-hole pair, is the sum of the maximum and 
the minimum values of these bands. (Remember that the Hartree-Fock eigenvalues at z/ = 1 
occur in opposite energy pairs.) It is the charge gap which is measured in the transport 
experiments. We see in Fig. 4 that, within the Hartree-Fock approximation, the effect of 
interactions on the charge gap is enormous. The hopping energy is here of order 0.01 (e^/£) 
while the band energy is one order of magnitude higher! However, we find that the Hartree- 
Fock charge gap is only weakly dependent on the parallel component of the magnetic field, 
in complete disagreement with experiment. This behavior is actually expected. The lowest- 
energy charged excitations are believed BJ6|j9|,p^ to be spin-textures and not the simple 
single-particle excitations whose energies are obtained from the Hartree-Fock calculations. 
The spin textures in double-layer systems are analogous to the skyrmion spin-textures which 
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are known to be the lowest energy charged excitations [|T7| , p!8[ of a single-layer system which 
is not spin-polarized. 

Fig. 5, shows the commensurate energy for various values of the bare hopping t. If we 
approximate the critical wave vector for the commensurate-soliton transition by that for the 
C -^ I transition (this is not a bad approximation as we saw above), we find the curve given 
in Fig. 6. Notice that because of the renormalization of t to t, there is always a critical 
value of Q at which this transition occurs. The angle that the magnetic field makes with 
the DQWS at the transition is given by tan(6'c) = QC/ {d/ 1) . In Fig. 6, we plot the critical 
angle for the transition at d/l =1. It is clear that the critical angle has a square-root 
dependence on t at small t. The same result was p[ obtained in the gradient approximation. 
Experimentally the critical tilt angle is identified with a feature in the dependence of the 
charge gap on tilt angle and it appears that the tilt angle has a linear dependence on t. 
Moreover, for given t the critical angle appears to be overestimated by both the Hartree- 
Fock approximation and by the gradient approximation when the stiffness coefficient, ps. 



see Eq.( |A14| )) is calculated in the Hartree-Fock approximation. |Q] For example, the critical 



angle that would be given for the parameters d = 1.877, t = 0.005 (see Fig.l (b) above) 
which are similar to those for the sample studied in Ref. |^, is approximately of 21° . This 
value is about 2.6 times greater than that of the experiment (8°). As discussed elsewhere [^ 
we believe that this quantitative discrepancy is due to quantum fluctuations not included 
in the Hartree-Fock approximation. For larger values of t these fluctuations are suppressed 
and we expect, in agreement with experiment, that critical tilt angles will agree more closely 
with the Hartree-Fock approximation. 

IV. RESPONSE FUNCTIONS OF THE COMMENSURATE AND SOLITON 

PHASES IN THE GRPA 

We now consider the calculation of the response functions of the DQWS in the presence 
of the tilted magnetic field. The total response is characterized in general by sixteen coupled 
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susceptibilities: 

Xtjkiiq, q', t) = -g < Tpij{q, r)pM(-q',0) >, (34) 

where pij = p— < pij >. It turns out (as in the ground-state calculation), that the equations 
of motion for the susceptibilities have a simpler form if we work with the shifted operators 
Prl,Plr instead of the unshifted ones. We define the 4x4 matrix 

r(q,q',r) = 

XRRRRifi, q', t) XRRLRi^l, q'+Q, t) XRRRLi<l, q' - Q, t) XRRLl{.^, q', r) 

x/?Li?,R(q + Q,q',^) x/?LLij(q + Q,q+Q,7-) XijL/?L(q+Q,q'-Q,^) XijLLL(q + Q,q',^) 

XLi?i?i?(q-Q,q',^) XL/?Li?(q - Q, q'+Q, 7-) XLi?i?L(q-Q,q'-Q,r) XLi?LL(q-Q,q',^) 

XLLiji?(q,q',T") XLLLi?(q,q' + Q,^) XLLi?L(q,q'-Q,^) XLLLL(q,q',^) 

(35) 



where again Q =Q^. 

The generalized random-phase approximation (GRPA) for the 2DEG in a strong mag- 
netic field with broken translational symmetry is discussed in detail in Ref. |]12|. The ap- 



plication of this formalism to the present problem presents no major difficulty although the 
algebra is knotty. Below we briefly sketch the derivation of the expressions we use for the 
numerical evaluation of the GRPA response functions. 
We first define the functions 

f{qx) = tSg^fi + {j^) Vdiqcc + Q) < pRLiqx) >, 

i^(q-q') = f{q.-q>-'^''-''^-''^^'^''^'5^,,,'^, 

F^(q-q') = F(q'-q), (36) 

K{q - q') = < PRL{q. - I) > e-fe-^^+Q)^^^V2 ^^^^^, ^ 

i^^(q-q') = ir(q'-q), 
and the 4x4 matrices 
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P{CL, q') 



^(q, q) 



^(q,q') 






-F^(q - q') 


i^(q-q') 





-F*(q-q') 








F(q-q') 


Ft(q-q') 








-i^^(q-q') 





Ft(q-q') 


-F*(q-q') 








-K^(q-q') 


^(q - q') 





-^*(q-q') 








K{q - q') 


irt(q_q') 








-ir^(q-q' 





irt(q-qO 


-ir*(q-q') 





H(q)<5q,q' 











1 


^^d(q + Q)5q,q' 














Vd(q-Q)^q,q 














H(q)'5q,q' 



i7(q,q') 



(37) 



(38) 



(39) 



(40) 



and 

K(q)<5q,q' K(q)5q,q' 





^ K(q)<^q,q' K(q)5q,q' 

Using h^{- ■ ■) = [H — jjN,- ■ ■] where H is the Hartree-Fock Haniihonian (Eq.(|^)), we 
find that the bare susceptibihty r° (corresponding to a one-loop approximation with the 
propagators evaluated in the HFA) satisfies the matrix equation 



Y: [1(00 + z6) 5q,q. - P(q, q")] r°(q", q', .;) = 5(q, q') 



(41) 



Here / is a unit matrix. We note that this bare susceptibility is independent of the value of 
qy since the propagators do not involve qy. This is no longer the case when we add bubbles 
and ladder diagrams. The summation of the ladder diagrams is given by the equation (F is 
the irreducible polarization) 



f (q, q', uj) = r°(q, q', u;) - E r°(q, q", a;)X(q", q'")? (q'", q', 



UJ 



(42) 



so that the full GRPA susceptibility is a sum of bubbles diagrams containing the irreducible 
polarization and is given by 



r(q,q',.;) = f(q,q',^)+ ^ f(q, q",o;)i7(q", q'")r(q'", q',^). 



(43) 



q",q" 



The Qy dependence enters only through the interaction H and X and so it is obvious that 

r(q, q', ^) = r(9x,gi; qy, ^)'^<?„<?,'- 

Because of the periodicity of the system, we must have r(q, q', u) = r(k+riQ, k+mQ, u), 
where n,m = 0, ±1, ±2, . . . and k is a vector in the first Brillouin zone of the modulated 



structure i.e. kx G 



— TT _7T_ 



and ky is unrestricted (the Brillouin zone is effectively one- 
dimensional). This implies that all elements in the above matrices are themselves matrices 
with respect to the indices n, m. Eqs. (^TD-(^3D can then be written in a purely matrix 
form and the structure of the GRPA becomes much more transparent. One must however 
be very careful not to interchange the order of these matrices since they do not commute! 
Combining these three equations, we arrive at 



[liu + iS) -R]T = S, 



(44) 



where 



R = P-S{X -H), 



(45) 



or, more explicitly 



R 



~F^ + K^X22 F-KX33 



F* + K*Xn - K*Hn + KH^ 








--^21 


^t - K^Xu + K^Hu - K'^Hu 








-^31 





— R^r> 


— -Rto 






(46) 



Looking at the form of the matrix R, it is quite clear that Eq.(^) can be further simplified. 
After some algebra, we find that 
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I{uj + i5f -A 



r = s, 



(47) 



where 



with 



and 



A 



All Ai4 

A22 A23 

A32 A33 

A*4 Ati 





All 


= R12R21 + -R13-R31, 






A22 


= R21R12 + -^21-^12' 






A33 


= -R31-R13 + -R31-R13, 






Al4 


= — -R12-R2I ~ -R13-R3I; 






A23 


= R21R13 + -R2I-R13) 






A32 


= i?3ii?i2 + -R31-R12' 




R12K* - RuK^ 


-{u 


+ i6)K^ {u + i6)K 


-R12K + R^K^ 


-{uj + i6)K* 


R21K 


^ + R*^^K^ -R2iK-R*^K* 


{uj + i6)K 


{u + i5)K^ 


RuK 


'' + Rl^K^ -R^iK-Rl.K* 


-{uj + t6)K^ 


RI2K* + Rl^K^ 


{u 


+ t6)K^ -{uj + i6)K* 


Rl^K - R\^K^ 



(50) 



In order to make connection with the pseudospin language used in Ref . [§] , we define the 
density and pseudospin operators by (remember that real spins are assumed fully polarized) 



(48) 



(49) 



(51) 



n = Prr + Pll, 

Sz = 2 (PRR " Pll) ■ 
For example, we can now construct the density response function as well as the longitudinal 
and transverse pseudospin response functions 



Xnn = XRRRR + XRRLL + XLLRR + XLLLL, 

XSzSz = 4 {Xrrrr " Xrrll — Xllrr + Xllll) 

X+ - = Xrllr- 



(52) 
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For the functions of Eq.(^), we need to solve for 



XL 



and 



XT 



/CRRRFl /xIxRLL 

Xllrr Xllll 

Xrllr Xrlrl 
Xlrlr Xlrrl 



(53) 



(54) 



With the definitions 



Ra 



Rto R 



R 



B 



-^12 



-R2I 
-R3I 



13 



-^13 



(55) 



-R 



21 



-R 



31 



(56) 



and 



K, 



K* -K 



(57) 



we arrive at the form of the equations which we have solved numerically: 



(58) 



[I{u + i5f-RARB\XL = RaKa 

[I{uj + i5f-RBRA]XT = RbK\. 
Notice that R12 is not simply related to R21 so that there is no simple relationship between 
Ra and Rb- These two matrices do not commute in the general case. We stress once again 
that Ra and Rb are not simple 2x2 matrices since each of their elements is itself an infinite 
matrix. Inverting Eq. (pSf), we easily arrive at 

XL = RA[I{UJ + I6f-RBRAV^KA, 

(59) 
XT = \I{Lo + i5f-RBRA] 'RbKI 

Where now the order of the matrices Ra and Rb are the same in both equations and the 
matrix RbRa is real. This shows clearly that all response functions do share the same poles 
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although each pole will in general have different residues for different response functions. 
The collective modes associated with the longitudinal response functions are thus the same 
as those associated with the transverse response functions. We can of course solve, in the 
same way, for all other response functions in Eg. (|35|) . 
If we diagonalize the matrix RbRa by 

(RbRa) U = Un, (60) 

where f2 is a diagonal matrix of the eigenvalues and U the matrix of eigenvectors, we can 
write 

When Qy = 0, all matrices are real and the above equations simplify considerably. In this 
case, we can easily show that 

^"'^ ^ °' (62) 

A. Sz Sz /A. RRRR ' 

There is thus no density response to potentials which are independent of the coordinate 
along the parallel magnetic field. 

The numerical procedure to compute the response functions and find the collective ex- 
citations is very simple. Once the order parameters are known from the HFA described in 
the previous section, we can compute the matrices Rb,Ra, and Ka- We then diagonalize 
the real matrix RbRa to find the eigenvectors U and eigenvalues Q. We can then compute 
the susceptibility from Eq.(pTD or, by following the pole with large weight as the momentum 
k varies in the Brillouin zone, find the dispersion relation of the collective modes. Before 
describing our numerical results for the soliton-lattice state, it is instructive to look at the 
response functions in the commensurate phase where the GRPA equations can be solved 
exactly. 
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A. Response functions in the commensurate phase 

In the commensurate phase, the only non-zero order parameter is < Phl(O) >= 1/2. It 
follows then that 

f{qx) = tRSq^fl, 

F(q,q') = tRe-'^5q,q., 
F^(q,q') = F(q,q'), 
^(q,q') = ie-'^<5q,q', 

ir^(q,q') = ir(q,q'), 



where we have defined 



and 



tR = t+-Vd{Q), 



^ 



(63) 



(64) 



It is very simple to solve for r(q, q',C(j) given in Eq.(^). First, we use the pseudospin 
transformation of Eq. (|5lD to define a new susceptibility matrix by 
rp(q,q',w) = 



/ 



Xn5.(q,q,w) 



Xnn(q,q',u;) Xn-(q,q'+Q,^) Xn+(q,q'-Q,^) 

x+n(q + Q,q',^) x+-(q+Q,q'+Q,^) x++(q + Q,q'-Q,^) x+5.(q + Q,q,^) 

x-n(q-Q,q',tu) X— (q-Q,q'+Q,tu) x-+(q-Q,q'-Q,^) x-5.(q-Q,q',^) 

X5.n(q,q',^) x5.-(q,q' + Q,c^) X5.+(q,q'-Q,^) X5.5.(q,q',^) j 

After some simple algebra, we find 

/ 

—2b{l~cos{2ip)) iujsm{ip) —iujsm{ip) ibsm{2(p) 

—iusm{ip) —a a —^cos{ip) 



V 



rp(q,q',c^) 



5q,q' 



{oo + z6y ~ ui{q) ^ 



iuj sin {if) 
-ifesin {2ip) 



a —a 

cos {if) ^ COS (ip) 



fcos((/.) 
;i + cos(2(/?)) 
(65) 



/ 
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where we have defined 

a = tfl + I (K(q) - H(q) - K(q) cos {2ip)) , 

h = tR-\ (\/,(q + Q) + VM - Q)) , (66) 

cc'g(q) = Aah. 

From Eg. (|55D , we see that the modulation of the pseudospin structure {i.e., of 0{X) ) 
due to the parallel magnetic field induces a coupling between the charge and spin response 
functions. In fact, as discussed in Ref. 0, this coupling is due to the magnetic field and 
technically arises because we have restricted the Hilbert space to the first Landau level 
only. In the restricted Hilbert space, pseudospin and charge operators no longer commute 
and it follows that modulations of the pseudospin texture are in general accompanied by 
modulations of the charge density. The presence of the parallel magnetic field is crucial, 
however, for this coupling to show up. In fact, the density response function is non-zero only 
in the presence of a parallel magnetic field. Without this field, the ground state consists of a 
fully filled Landau level of symmetric states. Since a scalar field coupling to the density does 
not change the pseudospin state, there is no way to excite an electron from the symmetric 
to the antisymmetric state. When the parallel field is added, however, symmetric and 
antisymmetric states are mixed and density wave excitations becomes possible. By density 
wave excitation, we mean here a collective excitation i. e. a pole of the two-particle Green's 
function. 

The coupling between pseudospin and density modes is interesting because the collective 
excitations of the system must now be considered as mixed pseudospin and density excita- 
tions. The coupling should also be important experimentally since external electromagnetic 
fields will tend to couple only to the total density of the two layers; in the presence of a 
parallel component of the magnetic field the density-density response function has a pole at 
the frequency of the pseudospin collective mode. In the soliton-lattice phase to which we 
turn next and, which has gap less collective modes, this coupling is also present. A similar 
situation is encountered in spiral magnets where the non-collinearity of the spin structure 
gives rise to a spin-charge coupling [jl9|. 
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B. Numerical results for the collective excitations 

In the commensurate phase, the various response functions x ^i^re very simple since they 
have only one real pole (no broadening). The pole of the bare response functions x° corre- 
sponds to the excitation of an electron (a particle-hole excitation) from the E^ to Ej^ band 
(Eq.p^D). When vertex corrections are introduced in the GRPA, this pole is renormalized 
into a dispersive collective excitation (a spin wave in the pseudospin language) and appears 
as a pole of x ■ These vertex corrections are important. To give an example, the gap calcu- 
lated from ojQ^q = 0) (Eq.(|5BD), at d/i = 1.0 and t/ (e^/eo^) = 0.01, is reduced by a factor 
six by the vertex corrections. Fig. 7 shows the dispersion relation of this spin wave, uJo{<i), 
for different values of the parallel magnetic field along the direction g^ = (the dispersion 
relation is not isotropic since the parallel magnetic field is set along the y axis). A soft mode 
develops as the parallel magnetic field is increased, but the value of Q at which the collec- 
tive mode of the commensurate state first goes soft is larger than the value of Q at which 
the commensurate to soliton-lattice phase transition occurs. The commensurate to soliton 
lattice phase transition is first order. No corresponding partial mode softening appears as a 
precursor of the commensurate to soliton lattice phase transition along the direction q^ = 0. 

The response functions of the soliton-lattice state are more interesting. Fig. 8 shows 
the imaginary part of the transverse response function x_|__(k + Q,k + Q,ct;) (a) in the 
HFA and (b) in the GRPA. The renormalization of the particle-hole excitation is more 
dramatic than in the commensurate phase. In addition to a Goldstone mode which appears 
at low energies, a number of additional sharp peaks suggestive of further collective modes 
appear in the response functions. For wave vectors in the x direction and parallel fields 
close to that of the commensurate to soliton lattice phase transition, the Goldstone modes 
are best thought of as compressional waves of the soliton lattice. For wave vectors in the 
y direction the Goldstone modes in this regime are best thought of as oscillations in the 
position of a particular phase slip as a function of the y coordinate. This Goldstone mode 
is expected because of the translational invariance of the underlying Hamiltonian which is 
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broken in the soliton-lattice state. The dispersion relation of the Goldstone modes is highly 
anisotropic as shown in Fig. 9 where it is computed for different values of Q along the 
two perpendicular directions g^ = and qy = 0. Close to the phase transition between 
commensurate and soliton-lattice states the solitons are very widely spaced and we should 
expect that the velocity of the Goldstone mode in the Qy = direction should be very small. 
In this direction (the direction perpendicular to -B||), the Goldstone mode is periodic with 
the periodicity of the Brillouin zone i.e. 27f/Ls. 

As we mentioned above, in the absence of tunneling, the ground state of the system is 
independent of parallel field, has spontaneous inter-layer phase coherence, and has Goldstone 
modes corresponding to this broken symmetry. The incommensurate state is the Hartree- 
Fock approximation for the broken symmetry state. For non-zero values of the tunneling 
amplitude in a parallel field, either the commensurate or the soliton-lattice state has a lower 
energy and the phase transition between these states occurs at a critical parallel field. For 
very strong magnetic fields, however, the soliton-lattice state asymptotically approaches the 
incommensurate state. The Goldstone modes of the soliton-lattice state, which are most 
naturally associated with broken translational symmetry, must crossover to the Goldstone 
modes of the spontaneous-phase-coherence state. These are easily obtained by solving GRPA 
equations at t = 0. (The calculations are very similar to those for the commensurate state). 
We find immediately that the response functions have a single pole at a frequency given by 

iOsM = ^[Va{q) - Vt{q) - V,{q) + Ki(0)][\/,(0) - VM]- (67) 

In Fig. 10 we compare the collective modes of the soliton-lattice state in the strong-field 
regime (where Qs ~ Q) with the collective modes of the spontaneously phase-coherent state 
in the absence of tunneling. As expected the collective modes are nearly identical. For 
the soliton-lattice state many collective modes appear at each wavevector; for {k^, ky) the 
energies are close to ujsc{\{kx + nQ, ky)\) where n is an integer. (The higher energy collective 
modes in Fig. 8 correspond to n ^ 0.) The broken translational symmetry of the soliton 
lattice state results in a folding of the collective modes of the t = ground state. At the 
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parallel field of Fig. 10 mode-mode coupling effects due to the interlayer hopping are weak; 
as the parallel field weakens and the phase transition is approached, the coupling becomes 
stronger and eventually the modes become more similar to those of the incommensurate 
state, except at wavelengths approaching the lengthening soliton-lattice period. 

At strong parallel fields, the effect of the tunneling on the spontaneously-coherent ground 
state is very much like the effect of a grating with wavevector Q. We are presently inves- 
tigating the possibility |^l| of using tunneling in a parallel field as a tunable grating which 



allows light to couple to finite wave vector modes of the spontaneously phase-coherent state 
in much the same way that artificially generated metallic gratings are used to couple to 
finite wave vector plasmons of two-dimensional electron systems. [^ 

Finally, we have verified that the pseudospin- charge coupling that we discussed in the 
commensurate phase also appears in the soliton phase. However, the pseudospin-charge 
coupling is very small. Typically, the density response function is more than one-thousand 
times smaller than the pseudospin response functions xzz or X+-- The weight of the Gold- 
stone mode in the density response function (as in all other response functions) vanishes as 
the wave vector of the excitation vanishes. Further studies will be needed to understand 
how this coupling affects the ground-state conductivity of the soliton lattice. 

V. CONCLUSION 

We have studied the commensurate-incommensurate phase transition which appears in 
a DQWS in a strong magnetic field when the field is tilted away from the normal to the 
planes of the two-dimensional electron gas in each well. Working in the Hartree-Fock ap- 
proximation, we have calculated the energy of the soliton-lattice state for different values of 
the tunneling parameter and parallel magnetic field and found that the soliton-lattice state 
is the ground-state of the system for all angles greater than the critical tilt angle. The dif- 
ference in energy between the soliton-lattice state and the incommensurate state is however 
very small. Close to the commensurate-soliton lattice transition, we find that the gradient 
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approximation described in Ref. P] is a good approximation for the soliton-lattice energy 
and for the modulation of the interlayer coherence phase. This approximation soon breaks 
down, however, as the parallel magnetic field is increased. 

Using the GRPA, we have also studied the collective modes of the soliton lattice. As 
expected from the broken translational symmetry of the soliton lattice, there is a Gold- 
stone mode branch which represents the phonons of the soliton lattice. Because of the 
one-dimensional character of the Brillouin zone, the Goldstone mode dispersion relation is 
periodic only in the direction perpendicular to the parallel magnetic field. Beside the Gold- 
stone mode, there are also higher-energy collective excitations in the response functions 
calculated in the GRPA. In the limit of strong parallel magnetic fields, we were able to 
relate these higher-energy excitations to the collective modes of the spontaneously- coherent 
state in the absence of tunneling at wavevectors displaced by multiples of Qx. In this regime, 
the effect of tunneling in a parallel field resembles that of a finite- wave- vector grating which 
mixes the different modes of the spontaneously coherent state. 

Finally, we have also shown that, except in the direction perpendicular to the parallel 
magnetic field, the collective modes in the system are in fact mixed pseudospin-charge exci- 
tations. Since all response functions are coupled in the soliton lattice, it follows that the den- 
sity response function has poles at exactly the frequency of these mixed pseudospin-charge 
excitations (including the Goldstone mode). More work is needed, however, to understand 
how this coupling could affect the conductivity of the soliton-lattice state. 
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FIGURES 
FIG. 1. Energy of the commensurate, incommensurate and soliton-lattice states as a function 

oiQl for (a) d/i = 1.0, t/ (eVeo^) = 0.01 and (b) d/l = 1.877, t/ (eVeo^) = 0.005. The dotted line 

corresponds to the energy of the commensurate phase calculated in the gradient approximation. 

FIG. 2. Behavior of 9{X) over one period of the soliton lattice in the HFA. (a) 
d/l = 1.0, t/(eVeo^) = 0.01,(5^ = 0.63 (just after the C ^ S" transition); (b) 
d/l = 1.0,t/(e2/eo^) = 0.01,(5^ = 1.5 (in the high parallel magnetic field limit). The dotted 
lines represent the sine-Gordon solution for 0[X) 

FIG. 3. Behavior of the ratio Qs/Q (soliton wave vector to parallel magnetic field wave vector) 
as a function of Qi for d/i = 1.0, t/ (e^/eo^) = 0.01. The line at Qs/Q = 1 is a guide to the eyes. 

FIG. 4. Shape of the upper band E^{X) over one period of the soliton lattice in the HFA for 
(a) d/i = 1.0, t/ {e^/eoi) = 0.01, Qi = 0.63 and (b) d/i = 1.0, t/ {e^/eoi) = 0.01, Qi = 1.5. 

FIG. 5. Energy of the commensurate state as a function of Qi for different values of the bare 
hopping parameter t. (1) t/ {e^/e^i) = 0.1; (2) t/ {e'^/e^i) = 0.01; (3) t/ {e'^/eoi) = 0.001. 

FIG. 6. Critical tilt angle at which the transition from commensurate to incommensurate states 
appears as a function of the bare hopping for d/i = 1.0. The full line is a fit to a square root 
dependence. 

FIG. 7. Dispersion relation of the spin wave mode in the commensurate phase for 
d/i = 1.0, t/ (e^/eo^) = 0.01 and different values of the parallel magnetic field. 

FIG. 8. Imaginary part of x+-(k + Q,k + Q,u;) for d/i = 1.0,t/ {e^/eoi) = 0.01, Qi = 0.72 
and wave vectors k^ = (0.1,0), (0.3,0). (a) HFA and (b) GRPA. The Goldstone modes appear as 
low-energy peaks in (b). 
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FIG. 9. Dispersion relations of the Goldstone mode at (a) dji = 1.0, t/ (e^/eo^) = 0.01, 
Qi = 0.72, and (b) djl = 1.0, tj {e^/eoi) = 0.01, Ql = 1.5 along the directions /c^ = and 
ky = 0. The dispersion relation is highly anisotropic, being in fact periodic only along ky = 0. 

FIG. 10. Comparison between the collective modes of the soliton lattice (•) at, 
t/ (e^/eo^) = 0.01, in the high-parallel magnetic field regime and the collective modes of the 
spontaneously-coherent state in the absence of tunneling (lines). The parameters are d/i = 1.0, 
QH. = 1.5 and ky = 0. For the soliton lattice, only the strongest six modes in the response function 

Xh are shown. For the spontaneously-coherent state, uJscikx + nQ) is plotted for n = —2, —1, ..., 4. 

The dotted line represents tUscikx)- 
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APPENDIX: THE GRADIENT APPROXIMATION 

In this appendix, we derive the Hartree-Fock approximation equation of motion for the 
single-particle Green's functions Gij{X, r) = — < Tcj^x(t)c]x(0) > instead of for its Fourier 
transform as in Eq.(|^). We then introduce the so-called gradient approximation developed 
in Ref. [Q. This approximation is very helpful in understanding how the system evolves 
from the commensurate to incommensurate phases in the presence of a parallel magnetic 
field. 

In X space, the Hartree-Fock Hamiltonian of Eq. (|^) is 

H = -tEx (e-^'^^'pRLiX) + e'Q^pLRiX)) 

(Al) 
- I {$) Ex,x' K{X - X') [< pn^iX) > Plr{X')+ < Plr{X) > Prl{X')] , 

where we again assume that the DQWS is unpolarized and neglect the constant first term 
in Eq.d). We have also defined the operator 

Pi AX) = clxCj,x, (A2) 

and the function 

K{X - X') = E Mq^) cos [qAX - X')] . (A3) 

The equations of motion for the single-particle Green's functions Gij{X,T) are 



{zuJn + p)GnR{X,Un) + S*{X)GLR{X,u;r,) = 1, 

{iuJn + P)Glr{X, UJn) + S{X)Grr{X, CJ„) = 0, 



(A4) 



where 



J{X) = i (^^^ Y.K{X-X') < pn,{X') >, (A5) 



and 



S{X)=te''^^ + J{X). (A6) 

Uncoupling Eqs.(|A^), we have 
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1 1 

2 I 2 



ILUn + fX - \S {X)\ lUJn + fJ. + \S {X)\ 



and 

1 / S{X) 



GLR{X,UJn) 



2\\S{X)\^ 



+ 



(A8) 



_iujn + jJ- - \S {X)\ iujn + jJ^+lS {X)\_ 

There are two bands with energies ± [^(X)!. At T = OK, one of these bands must be 
entirely filled with the other one empty in order to insure that the filling factor of the 
DQWS be exactly one. Thus, 

or, 

This is the Hartree-Fock equation for the order parameter in X space. We have devised 
an iterative scheme to solve this equation and at the same time minimize the Hartree-Fock 
ground-state energy per particle which is given by 

e = --y: cos [e{x) - QX] - T^(4) E K{x - X') COS {e{x) - e{x')). (aii) 

We obtain results which are identical to those obtained using the approach described in the 
body of this paper. 

It is easy to see that Eq.( [Al(J| ) can be obtained by minimizing the energy with respect 
to the phase i.e. with jg^ = 0. The Hartree-Fock solution is thus an extremum of the 
Hartree-Fock energy as expected. 

The function K{X — X') is a rapidly decreasing function of \X — X'\. If the angle 0{X) 
does not vary too rapidly in space, this suggests that we can expand the cosine function in 
Eq.( [All| ). We get in this way 

^ - -i (t) ^^(°' - 5 S COB m) - QAI + 5^ E (^) ^ (A12) 
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where 7 = (77) YjX K{X)X'^. Going to the continuum hmit with YjX ■ ■ ■ — ^ 2^ -^ '^■^ • • •' 
we have finally for the total energy of the system: 



E^c:,-U^]VM+Ly I dX 



where e{X) = 9{X) - QX and 



1 fd9{X) 



:Ps 



\ dX 



Q 



27r£2 



cos(^(X)) 



(A13) 



(A14) 



Ps= {ii)TE^^IdXK{X)X^^ 

\eoiJ 87r£2 dQ2 q^q ■ 

If we now minimizes this total energy with respect to 0{X), we find the well-known 
sine- Gordon equation 

d'^e{X) t 



dX^ 27ipsi 
which admits the kink soliton as a solution: 

e{X) =4tan- 



sin(e(X)), 



(A15) 



g V '^^ps^ 



tX 



(A16) 



If we insert Eq.( [A16D in Eq.( [A13| ), we find that the energy per particle in the presence of a 
single soliton is 



Ess ^Ec + Lyps 



\ 2npsi^ 



2'kQ 



(A17) 



The first term in the right-hand-side of this equation is the energy per particle in the 
commensurate phase given in Eq.(^) while the second term is the change in this energy due 
to the presence of the soliton. When this second term becomes less than zero, we expect a 
phase transition in the system to occur since then the energy to create one soliton becomes 
negative. This happens at a parallel magnetic field given by 



B\ 



{c^s) 



Bi 



2/ 

nd 



2t 



\ -^Ps' 



(A18) 



If we define the width of the soliton as 
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e = ^^, (A19) 



2-KpsP 

then, just at the transition, we have 



ic = -^. (A20) 



The critical magnetic field given in Eq.( |A18 ) should be compared with the critical field 



for the transition from the commensurate to the incommensurate phase. In the gradient 
approximation, this critical field is given by given by 



Finally, for a periodic ground state with regularly spaced solitons, the number of solitons 
is given by Ns = L^/Ls where Ls is the distance between two solitons. If we neglect the 
soliton-soliton interaction, then the energy of the configuration is approximately given by 

Ec^Ec + PsiQsi) m - 2vrQ^] , (A21) 

where Qs = 2n/Ls. For Q > 4/(7r^) this energy is always lowered by increasing Qs and 
hence the density of solitons. Repulsive interactions between the solitons must be included 
to determine the equilibrium value of Qs- 
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